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We study transport properties and also local density of states in a clean superconductor-normal 
metal-superconductor Josephson junction with triplet /—wave superconductors based on the Eilen- 
berger equation. Effects of thickness of normal metal and misorientation between gap vector of 
superconductors on the spin and charge currents are investigated. Spin current, which stems from 
the misalignment of d-vectors, in the absence of charge current for some values of phase difference 
between superconductors is found. We also find unconventional behavior of the spin current associ- 
ated with 0-7r transition. The misalignment of d-vectors also gives rise to a zero energy peak in the 
density of states in the normal metal. 
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I. INTRODUCTION 



and density of states in a triplet SNS junction. 



Superconductivity with spin-triplet pairing symme- 
try which has been observed in a series of experiments 
in Sr 2 RuOj™ UPt 3 and some other heavy fermion 
complexes^^ attracts a lot of interests in the condensed 
matter physics. The spin triplet /—wave pairing sym- 
metry has been considered for unconventional supercon- 
ductivity in UPt3 heavy-fermiorP^HlS]. Superconductiv- 
ity has been observed in UPt3 compound in two impor- 
tant states, A— phase and B— phase. Although different 
phases correspond to the different symmetries of the su- 
perconducting gap on the Fermi surface, all of them be- 
long to spin triplet. The /—wave order parameter has a 
more complex dependence on azimuthal angles in com- 
parison with the well-known p— wave order parameter for 
superfluid phases of 3 He. The Fermi surface of UPt3 is 
highly complex with five or six complicated shape sheets. 
While the pairing mechanism of UPt3 is unknown, a large 
number of experimental and theoretical works have been 
done on its thermodynamic and transport property. The 
Josephson junction between triplet superconductors is an 
accurate method to study physics in triplet superconduc- 
tors. In superconductor-normal metal-superconductor 
(SNS) junction, the Josephson current flows owing to the 
proximity effect. 

Recently, proximity effect in unconventional supercon- 
ductor and normal metal junction has attracted much 
attentiorPSHUl xhe proximity effect in junctions between 
triplet superconductors and normal metals in the dif- 
fusive regime have been investigated theoretically^" 2 ^. 
It is found that the proximity effect is enhanced by the 
mid-gap Andreev bound state formed at the interface be- 
tween triplet superconductor and the normal metal due 
to the sign change of the pair potential.^ As a result, the 
density of states in the normal metal has a zero energy 
peak rather than gap. To investigate the proximity ef- 
fect, the quasiclassical Green's function method is quite 
useful. In this paper, we apply the quasiclassical Eilen- 
berger equation^ to calculate charge- and spin-currents 



The SNS Josephson effect is an interesting issue be- 
cause of the industrial applications such as supercon- 
ducting quantum interference devices (SQUID) 26 * 27 !. DC 
Josephson junction has a persistent current by gradient 
of the macroscopic phase of superconductivity 28 . While 
the charge current is more known property in Joseph- 
son junction, spin current is more interesting because of 
applications in spintronics. The effect of phase differ- 
ence on the spin current in the spin -tripl et SIS Josephson 
junction has been studied recentljE2E21 Authors of the 
papeil^l obtained a polarized dissipation-less supercur- 
rent of spins. It is shown that the current-phase depen- 
dencies are quite different from those of th e jun ction be- 
tween conventional s-wave superconductor s] 30 ! 31 ! and high 
temperature d-wave superconductorjP^. 

In the present paper, we investigate spin and charge 
Josephson currents, and the density of states in junction 
between misorientated crystal of triplet /—wave super- 
conductors which sandwich the mesoscopic normal metal 
layer. Triplet /-wave pairing symmetry is considered, 
for example, for superconductivity in UPt3 crystalJ^D. 



The organization of the rest of this paper is as follows. 
In Sec(TT] the quasiclassical equations for Green's func- 
tions are presented. Green's functions are obtained in 
Sec lIIII The obtained formulas for the Green's functions 
are used to calculate the charge and spin current densi- 
ties and also density of states in the normal region. In 



Sec. IV numerical results of charge and spin currents in 
the normal metal between superconductors are analyzed. 
Effects of normal layer thickness between superconduc- 
tors and misorientation between d-vectors of supercon- 
ductors on the charge and spin currents, and density of 
states have been investigated. The paper will be finished 
with the conclusions in SecfVl 
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FIG. 1: (Color online) Scheme of a flat interface between two 
superconducting bulks. Here, 1, 2 label the left and right 
half-spaces. Two superconducting bulks which are separated 
by a layer of normal metal have a misorientation as much as 
a. Two different geometries are corresponding to the differ- 
ent orientations of the crystals in the right and left sides of 
the interface. In geometry (i) the 6c— plane in the right side 
has been rotated around the a— axis by a. In geometry (ii), 
the ab— plane in the right side has been rotated around the 
c— axis by a. For both of geometries (i) and (ii), we consider 
a rotation only in the right superconductor and the crystal- 
lographic a— axis in the left half-space is selected parallel to 
the partition between the superconductors (a;— axis). 



II. FORMALISM AND BASIC EQUATIONS 

In this section, we consider a clean normal metal such 
as Copper between two misorientated /—wave supercon- 
ductors. A normal metal layer with the thickness of I is 
sandwiched by two triplet superconductors. Interfaces of 
normal metal and superconductors have been considered 
as totally transparent. For the case I > Ap we can use 
the ballistic quasiclassical Eilenberger equation^ 



hw F ■ V<? + e m a 3 +iA,g 



= 0, 



(1) 



and the normalization condition gg = I, where e m — 
TtkBT(2m + 1) are discrete Matsubara energies with m = 
0, 1,2..., T is the temperature, vp is the Fermi velocity 
and <7 3 = T 3 Cg)/ in which f 3 is the Pauli matrix in particle- 
hole space, dj (j = 1, 2, 3) denote Pauli matrices in spin 
space in the following. 

The Matsubara propagator g can be written in the 
standard form: 

- _ ( glff + 91 g2<™<T 2 \ , 2 s 

\ *02g3<T <X 2 gl 0-0-2 - 91 )' 

where the matrix structure of the off-diagonal self energy 
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FIG. 2: (Color online) The normal component (z component) 
of the char ge c urrent versus the phase difference <f> for the 
axial state (12 1, T = 0.05Tc, different thicknesses of normal 



metal and different misorientations between superconductors. 
Panel (a) is for geometry (i) and panel (b) for geometry (ii). 
Currents are given in units of jo,e = f eiV(0)«F A(0). 



A in the Nambu space is 



A = 





ia2<i*a 







(3) 



In this paper, we focus on the unitary states, (d x d* = 
0). Also, we use the Eilenberger equation for d = 
in the normal metal region (0 < z < I). From the 
Eilenberger equation it is clear that Green's functions 
in the normal metal, giN and gijsr, are constant for 
< z < I. Solutions of Eq.([T| has to satisfy the condi- 
tions for Green's functions at the bulk of superconductors 

V e m + I d 2,l| 

In addition, solutions of Eqs. satisfy the continuity 
conditions at the interfaces between metal and supercon- 
ductors (z = 0, 1) for all quasiparticle trajectories. 
Here, as in the paper^l, a simple step-like non-self- 
consistent model of the constant order parameter up to 
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FIG. 4: (Color online) The normal component of the spin 
current (j 3X )z versus the phase difference <f> for the ax- 
ial state ( |I2[ |, geometry (i) and the different thickness of 
metal and different misorientation, T = 0.05Tc and the 
different thicknesses of layer normal metal between super- 
conductors. Spin currents {j ay )z, (jsz)z for geometry (i) 
and (j ax )z,(jsy)z,(jsz)z for geometry (ii) are absent be- 
cause of Eq. Q. Currents are calculated in units of jo,s = 
?rhN (0)v F A(W). 



III. ANALYTICAL RESULTS OF GREEN'S 
FUNCTIONS 



FIG. 3: (Color online) The normal component (z component) 
of the charge cu rrent versus the phase difference <f> for the 
planar state (13 1, T = 0.05Tb, different thicknesses of layer 



normal metal between superconductors and different misori- 
entations between superconductors. Panel (a) is for geometry 
(i) and panel (b) for geometry (ii). 



the interfaces is considered: 



&{z,\f) = 



di(vF )e+'^/ 2 z < 
< z < I 



d 2 (v F )e- 



(4) 



> I 



where <j) is the external phase difference between the gap 
functions of superconducting bulks. We assume that the 
order parameter does not depend on coordinates and in 
each superconductors equals to its value far from the in- 
terface in the left or right bulks. For such a model, the 
current-phase dependence of a Josephson junction can 
be calculated for specific model of /—wave pairing sym- 
metry. We believe that under this assump tions our re- 
sults describe the real situation qualitativel y^ 4 * 35 !. In the 
framework of such model, the analytical expressions for 
the charge and spin currents and density of states can be 
obtained for an arbitrary form of the gap vector. 



The solution of Eilenberger equations allows us to cal- 
culate the charge and spin current densities in normal 
metal. The expression for the charge current is: 

j e (r) = 2meTN (0) £ (v F5l (v F , r, e m )) , (5) 

771 

that for the spin current is: 

j Si (r) = inhTN (0) ^ (v F (e ig i (v F , r, e m ))) (6) 

m 

and that for local density of states at the energy E is: 

N(E, r) = N (0) (Re 9l (v F , r, s m -> iE + 8)) (7) 

where (...) stands for averaging over the directions of 
electron momentum on the Fermi surface vp, N (0) is 
the electron density of states at the Fermi surface and 
e 4 = (x,y,z)- 

The calculated Green's functions in the normal metal 
are the following: 



9in 



,1N 



rj(A - B) 



A + B + 2\d 1 -d%\ 2 



27 ? dl-d 2 (dixd$) 
A + B + 2\d x -d*| 2 



(8) 



(9) 
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FIG. 5: (Color online) The normal component of the spin cur- 
rent (jsx)z versus the phase difference cf> for the planar state 
(13 1, geometry (i), T = 0.05Tc, different thicknesses of layer 
normal metal between superconductors and different misori- 
entations between superconductors. Note that spin currents 
(jsy)z and (jsz)z vanish. 
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FIG. 6: (Color online) The normal component of the spin 
current (j S z)z versus the phase difference 4> f° r the planar 
state ( |13p , geometry (ii), T = 0.057c , different thicknesses 
of layer normal metal between superconductors and different 
misorientations between superconductors. Spin currents of 
(jsx)z and (jsy)z are absent. 
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FIG. 7: (Color online) Critical currents of charge and spin 
(jsx)z as a function of misorientation angle for axial state 
Q, Z/£ = 0.5, T = 0.05T C . Solid lines are for geometry (ii) 
and dotted lines are for geometry (i). Here and in Fig[8] we 
define j ce = Max^,j e {<j>) = j e (<f>*) and j cs = j s {<t>*)- 
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FIG. 8: (Color online) Critical charge and spin c urren ts as a 
function of misorientation angle for planar state (13 1, £/£ = 
0.5, T — 0.05T C . Solid lines are for geometry (ii) and dotted 
lines are for geometry (i). 



Using these Green's function, we obtain the charge and 
spin current through the SNS Josephson junction and 
also the density of states of system in the normal region 
< z < I. 



where r\ = sgn (v x ), Q„ = J + |d„ 



+2e I 

A = d{- d 2 (e m + ??^i)(£ m + ^2) exp( "' ) (10) 

\Vz\ 



and 



— 2e m l 

B = di • d 2 (e m - r)£li)(£ m - 77^2) exp( ). (11) 

\v z \ 



IV. NUMERICAL RESULTS OF CURRENTS 
AND DENSITY OF STATES 

In this section we calculate charge and spin currents 
in SNS junctions having in mind UPt3 that is a /—wave 
superconductor. Two famous models for order parameter 
of UPt3 are considered. The first is the axial stated 



d(T,Vf) = A(T)z k z (k x 



iky) , 



(12) 
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FIG. 9: (Color online) Local density of states in the normal 
metal for a = 0, £/£ = 0.5, T = 0.05T C and for both symme- 
tries (axial and planar states) and different phases <j>. 
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FIG. 10: (Color online) Local density of states in the normal 
metal for both geometries with a — <j> = ^, T — 0.05T C , 
and different thicknesses of normal metal I. Panels (a) and (b) 
are for axial state ( |12[ | and panels (c) and (d) are for planar 
state (13 1. Left plots are for geometry (i) and right plots are 



f2k X ky] 



(13) 



for geometry (ii). 

and the other is the planar stated 

d(T,vj,) = A(T)fc a [x (k 2 x -k 2 y ) 

The coordinate axes x, y,z are chosen along the 
crystallographic axes a, b, c in the left side of Figjl] 
The function A (T) in Eqs.Q and ([T3j) describes 
the dependence of the order parameter d on the 
temperature T. Using Green's function in Eqs.([8]) 
and ([9]), we have calculated charge and spin currents, 
and the density of states in /—wave Josephson junction 



FIG. 11: (Color online) Local density of state in the normal 
metal for Z/£ = 0.5, T = 0.05T C , = § , and two different 
misorientations. Panels (a) and (b) are for axial state ( 12 1 
and panels (c) and (d) for planar state (13 1. Left plots are for 
geometry (i) and right plots are for geometry (ii). In panel 
(b), plots of a = j and a = ^ (which are coincident) are 
shown by a solid line. 



with pairing symmetry of Eqs. ( 12 1 and ( 13 1 numerically. 

Two superconducting bulks may have misorientation 
which has been created by geometries (i) and (ii) in Fig. 
[T] For two geometries and specific models of /—wave 
pairing symmetry we have plotted the currents in terms 
of phase difference in Figs(2][6j The obtained currents 
are periodic functions of phase difference between su- 
perconductors with a period of 2ir. It is clarified that 
by increasing the thickness of normal layer, amplitude 
of currents decreases, which is understandable because 
by increasing thickness of the normal metal, quantum 
coherency of macroscopic phases between left and right 
superconductors decreases. The current-phase relations 
are qualitatively different f rom th e case of s— wave^ and 
d— wave Josephson junction^ 33 ! 34 !. In particular, an inter- 
esting case is the zero of charge current at finite phase 
4>o in Fig{2] (b) which occurs for the most of junctions 
exactly at <j> = or <j) = ?r- On the other hand, for finite 
misorientation between left and right gap vectors even at 
the zero phase difference we obtain finite charge current. 
This is because according to Eq. ( 12 ) in geometry (ii) 
misorientation a plays role of phase difference -rotation 
by a gives the phase factor exp(2za). Also, in Figs. 4]|6 



it is found that the spin current of Josephson junction 
between two triplet superconductors has been created by 
misorientation between gap vectors of superconductors^. 
For the case of geometry (ii) and axial state (12 1, spin 
current is absent because both of gap vectors are in the 
same direction. Also, by using Eq. we find that 

spin currents (j sy ) z and (j S z)z f° r the axial state (12) 
and also the planar state (13) in geometry (i), and spin 



6 



currents (j S x)z an d (j sy )z for the planar state (13) in 
geometry (ii) are absent. In contrast to the charge cur- 
rent, the spin current has a large value at <j) = tt (or at 
= for large misorientation angle). It is seen that 
charge current are generally odd function but spin cur- 
rent are even function of the phase difference because 
charge current is odd in time reversal while spin current 
is even. Since spin current is an even function of phase, 
its derivative becomes odd in phase and hence should be 
zero at = and tt. Thus, the spin current has a lo- 
cal maximum or minimum at <j) — and tt. We also see 
that for some phases as <fr — tt the spin current exists but 
charge current disappears. This is a purely spin trans- 
port in the absence of charge current. This effect can not 
be observed in singlet superconducting junction like con- 
ventional or high T c superconductors. Only in the case 
of superconductor-ferromagnet-superconductor junction 
with inhomogeneous magnetization this spin current is 
presentPn 

The critical currents as a function of the misorienta- 
tion angle are shown in Figs[7] and [8] Critical charge 
and spin currents, defined as j ce = Max t j > j e {4>) = j e {4>*) 
and j cs = j s (<P*), respectively, are plotted as a function 
of misorientation angle a for specified temperature and 
normal layer thickness (T = &, I 
the misorientation, we see the — 
the anisotropy of the gapP^ Also, at the — tt transi- 
tion point for both symmetries in geometry (i) we have 
a jump of spin current due to the jump of the phase as 
found in RefP^but in geometry (ii) there aren't any — tt 
transition and hence jump of spin current. 

Finally, in Figs 9]flT we have plotted the density of 
states in the normal region. It is known that proximity 
effect in the SNS junction strongly changes the density 
of states in the normal metal. In our calculations, we ob- 
tained Andreev bound states in the junctions as shown in 
Figs|9]|TT) This is understandable, by interference effect 
of electrons and holes in the normal metal. An electron 
will be reflected by Andreev reflection as a hole. This 
process is also repeated on the other SN interface. Then, 
interference of holes and electrons increases the effective 
Andreev reflection probability and proximity effect at the 



= |). By changing 
tt transition due to 



low energies. 

In Figj9j we obtain a mini-gap like s— wave and d— wave 
superconductors^. Exactly at — tt we obtain a peak at 
the Fermi energy. Here, like Josephson junction between 
other types of unconventional superconductors^, we ob- 
tain the finite density of states for energies inside the gap 
(—A < E < A). The results of a finite misorientation 
and different I are shown in Fig. [TUJ The density of states 
gets smeared by increasing I. In Fig. [IT] we plotted the 
density of states for both of geometries and both of sym- 
metries. It is observed that by increasing the misorienta- 
tion two peaks are gradually merged into one and there 
will be a peak at zero energy at a = J . The increase of 
the density of states at low energy corresponds to that of 
the spin current. We find that effect of thickness of nor- 
mal layer is only quantitative but misorientation changes 
currents and the density of states qualitatively. 



V. CONCLUSIONS 

In this paper, we have investigated transport prop- 
erties and also local density of states in triplet SNS 
Josephson junction with two /—wave superconductors. 
A spin-polarized current normal to the interface is inves- 
tigated theoretically. It is observed that normal layer de- 
creases both of charge and spin currents. It is also found 
that misorientation of d- vector between left and right su- 
perconducting bulks produces the spin-polarized current 
which can flow even in the absence of charge current. 
We also unveiled a jump of the spin current associated 
with the 0-7T transition. In addition, density of states 
have been calculated for various phase difference, normal 
metal thickness and misorientation. Normal metal thick- 
ness has only quantitative effect on the mid-gap states. 
Misorientation of the <i-vectors drastically changes the 
junction property. 

A generalization of calculations of this paper to 
superconductor-ferromagnet-superconductor junctions is 
an interesting subject. We will do this generalization in 
the next work. 
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